Set up
suppressPackageStartupMessages({
library(tidyverse)
})
Read in data and process
pbta_df <- readr::read_tsv(pbta_file, guess_max = 100000, show_col_types = FALSE) %>%
select(Kids_First_Participant_ID, Kids_First_Biospecimen_ID, cg_multiple, cg_id, cgGFAC, tumor_descriptor)
genomic_paired_df <- readr::read_tsv(genomic_paired_file, guess_max = 100000, show_col_types = FALSE) %>%
left_join(pbta_df, by = c("Kids_First_Participant_ID"))
descriptors_df <- genomic_paired_df %>%
group_by(Kids_First_Participant_ID) %>%
summarize(descriptors = paste(sort(tumor_descriptor), collapse = ", "),)
# Vector to order timepoints
timepoints <- c("Diagnosis", "Progressive", "Recurrence", "Deceased", "Second Malignancy", "Unavailable")
# Create df to use for plots
df_plot <- genomic_paired_df %>%
left_join(descriptors_df, by = c("Kids_First_Participant_ID", "descriptors")) %>%
mutate(td_cgGFAC = case_when(grepl("Deceased", tumor_descriptor) ~ "xDeceased",
TRUE ~ tumor_descriptor),
log10_tmb = abs(log10(tmb)),
cg_id_kids = paste(cg_id, Kids_First_Participant_ID, sep = "_"),
cg_id_kids = str_replace(cg_id_kids, "/", "_"),
cg_id_kids = str_replace(cg_id_kids, "-", "_"),
cg_id_kids = str_replace_all(cg_id_kids, " ", "_"),
tumor_descriptor = factor(tumor_descriptor),
tumor_descriptor = fct_relevel(tumor_descriptor, timepoints))
# Let's count number of samples
count_df <- df_plot %>%
group_by(tumor_descriptor, cg_id, Kids_First_Biospecimen_ID, Variant_Classification) %>%
dplyr::count(cg_id)
Define parameters for plots
# Read color palette
palette_df <- readr::read_tsv(palette_file, guess_max = 100000, show_col_types = FALSE)
# Define and order palette
palette <- palette_df$hex_codes
names(palette) <- palette_df$color_names
Alterations per timepoint
# Define parameters for function
x_value <- count_df$tumor_descriptor
title <- paste("Variant types in PBTA cohort", sep = " ")
# Run function
fname <- paste0(plots_dir, "/", "Alteration_type_timepoints_barplots.pdf")
print(fname)
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/Alteration_type_timepoints_barplots.pdf"
p <- create_stacked_barplot_variant(count_df = count_df, x = x_value, palette = palette, title = title)
pdf(file = fname, width = 6, height = 6)
print(p)
dev.off()
quartz_off_screen
2

Alterations per timepoint in each cancer type
cg_id_id <- as.character(unique(count_df$cg_id))
cg_id_id <- sort(cg_id_id, decreasing = FALSE)
cg_id_id
[1] "Adamantinomatous Craniopharyngioma" "Atypical Teratoid Rhabdoid Tumor"
[3] "Chordoma" "Choroid plexus carcinoma"
[5] "CNS Embryonal tumor" "Craniopharyngioma"
[7] "Diffuse midline glioma" "Dysembryoplastic neuroepithelial tumor"
[9] "Embryonal tumor with multilayer rosettes" "Ependymoma"
[11] "Ewing sarcoma" "Ganglioglioma"
[13] "Glial-neuronal tumor" "Hemangioblastoma"
[15] "High-grade glioma" "Low-grade glioma"
[17] "Malignant peripheral nerve sheath tumor" "Medulloblastoma"
[19] "Meningioma" "Neuroblastoma"
[21] "Neurofibroma/Plexiform" "Pilocytic astrocytoma"
[23] "Rosai-Dorfman disease" "Schwannoma"
# Define parameters for function
x_value <- count_df$Kids_First_Biospecimen_ID
title <- paste("Variant types in PBTA cohort across cancer groups", sep = " ")
# Run function
fname <- paste0(plots_dir, "/", "Alteration_type_timepoints_cg_id_barplots.pdf")
print(fname)
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tmb-vaf-longitudinal/plots/Alteration_type_timepoints_cg_id_barplots.pdf"
p <- create_stacked_barplot_variant_cg_id(count_df = count_df, x = x_value, palette = palette, title = title)
pdf(file = fname, width = 30, height = 18)
print(p)
dev.off()
quartz_off_screen
2

Alterations per timepoint in each cancer type defined by cgGFAC
df_plot_cgGFAC <- df_plot %>%
filter(!is.na(timepoints_models)) %>%
arrange(timepoint_cgGFAC_n) %>%
group_by(tumor_descriptor, cgGFAC, timepoint_cgGFAC_n, Kids_First_Biospecimen_ID, Variant_Classification) %>%
dplyr::count(timepoint_cgGFAC_n)
cgGFAC_id <- as.character(unique(df_plot_cgGFAC$cgGFAC))
cgGFAC_id <- sort(cgGFAC_id, decreasing = FALSE)
cgGFAC_id
[1] "ATRT" "DMG" "HGG" "LGG" "Other"
# Define parameters for function
x_value <- df_plot_cgGFAC$Kids_First_Biospecimen_ID
title <- paste("Variant types in PBTA cohort across cgGFAC", sep = " ")
# Run function
p <- create_stacked_barplot_variant_cgGFAC(count_df = df_plot_cgGFAC, x = x_value, palette = palette, title = title)

Alterations per timepoint in each cancer type and timepoint
model
tm_df_plot <- df_plot %>%
filter(!is.na(timepoints_models)) %>%
group_by(tumor_descriptor, cg_id, timepoints_models, Kids_First_Biospecimen_ID, Variant_Classification) %>%
dplyr::count(timepoint_cgGFAC_n)
tm <- as.character(unique(tm_df_plot$timepoints_models))
tm <- sort(tm, decreasing = FALSE)
tm
[1] "Dx-Dec" "Dx-Pro" "Dx-Pro-Dec" "Dx-Pro-Rec" "Dx-Pro-Rec-Dec" "Dx-Rec" "Dx-Rec-Dec" "Dx-SM"
[9] "Pro-Dec" "Pro-Rec" "Pro-Rec-Dec" "Rec-Dec" "Rec-SM"
# Loop through variable
for (i in seq_along(tm)){
print(i)
df_sub <- tm_df_plot %>%
filter(timepoints_models == tm[i])
# Define parameters for function
x_value <- df_sub$Kids_First_Biospecimen_ID
title <- paste(tm[i])
# Run function
p <- create_stacked_barplot_variant_cg_id(count_df = df_sub, x = x_value, palette = palette, title = title)
}
[1] 1
[1] 2
[1] 3
[1] 4
[1] 5
[1] 6
[1] 7
[1] 8
[1] 9
[1] 10
[1] 11
[1] 12
[1] 13













sessionInfo()
R version 4.2.3 (2023-03-15)
Platform: aarch64-apple-darwin20 (64-bit)
Running under: macOS Ventura 13.5.2
Matrix products: default
LAPACK: /Library/Frameworks/R.framework/Versions/4.2-arm64/Resources/lib/libRlapack.dylib
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] grid stats graphics grDevices utils datasets methods base
other attached packages:
[1] ggthemes_4.2.4 lubridate_1.9.2 forcats_1.0.0 stringr_1.5.0 dplyr_1.1.2 purrr_1.0.1 readr_2.1.4 tidyr_1.3.0
[9] tibble_3.2.1 ggplot2_3.4.2 tidyverse_2.0.0
loaded via a namespace (and not attached):
[1] tidyselect_1.2.0 xfun_0.39 bslib_0.5.0 carData_3.0-5 colorspace_2.1-0 vctrs_0.6.3 generics_0.1.3
[8] htmltools_0.5.5 yaml_2.3.7 utf8_1.2.3 rlang_1.1.1 pillar_1.9.0 jquerylib_0.1.4 ggpubr_0.6.0
[15] glue_1.6.2 withr_2.5.0 bit64_4.0.5 lifecycle_1.0.3 munsell_0.5.0 ggsignif_0.6.4 gtable_0.3.3
[22] ragg_1.2.5 evaluate_0.21 labeling_0.4.2 knitr_1.43 tzdb_0.4.0 fastmap_1.1.1 parallel_4.2.3
[29] fansi_1.0.4 broom_1.0.5 scales_1.2.1 backports_1.4.1 cachem_1.0.8 vroom_1.6.3 jsonlite_1.8.7
[36] abind_1.4-5 systemfonts_1.0.4 farver_2.1.1 bit_4.0.5 textshaping_0.3.6 hms_1.1.3 digest_0.6.33
[43] stringi_1.7.12 rstatix_0.7.2 rprojroot_2.0.3 cli_3.6.1 tools_4.2.3 magrittr_2.0.3 sass_0.4.7
[50] crayon_1.5.2 car_3.1-2 pkgconfig_2.0.3 timechange_0.2.0 rmarkdown_2.23 rstudioapi_0.15.0 R6_2.5.1
[57] compiler_4.2.3
LS0tCnRpdGxlOiAiQ2xhc3NpZmljYXRpb24gb2YgVmFyaWFudHMgYWNyb3NzIHBhaXJlZCBsb25naXR1ZGluYWwgc2FtcGxlcyBpbiB0aGUgUEJUQSBDb2hvcnQiCmF1dGhvcjogJ0FudG9uaWEgQ2hyb25pIDxjaHJvbmlhQGNob3AuZWR1PiBmb3IgRDNCJwpkYXRlOiAiMjAyMyIKb3V0cHV0OgogIGh0bWxfbm90ZWJvb2s6CiAgICB0b2M6IFRSVUUKICAgIHRvY19mbG9hdDogVFJVRQotLS0KCiMgU2V0IHVwCmBgYHtyIGxvYWQtbGlicmFyeX0Kc3VwcHJlc3NQYWNrYWdlU3RhcnR1cE1lc3NhZ2VzKHsKICBsaWJyYXJ5KHRpZHl2ZXJzZSkKfSkKYGBgCgojIERpcmVjdG9yaWVzIGFuZCBGaWxlIElucHV0cy9PdXRwdXRzCmBgYHtyIHNldC1kaXItYW5kLWZpbGUtbmFtZXN9CiMgRGV0ZWN0IHRoZSAiLmdpdCIgZm9sZGVyIC0tIHRoaXMgd2lsbCBiZSBpbiB0aGUgcHJvamVjdCByb290IGRpcmVjdG9yeQojIFVzZSB0aGlzIGFzIHRoZSByb290IGRpcmVjdG9yeSB0byBlbnN1cmUgcHJvcGVyIHNvdXJjaW5nIG9mIGZ1bmN0aW9ucyAKIyBubyBtYXR0ZXIgd2hlcmUgdGhpcyBpcyBjYWxsZWQgZnJvbQpyb290X2RpciA8LSBycHJvanJvb3Q6OmZpbmRfcm9vdChycHJvanJvb3Q6Omhhc19kaXIoIi5naXQiKSkKYW5hbHlzaXNfZGlyIDwtIGZpbGUucGF0aChyb290X2RpciwgImFuYWx5c2VzIiwgInRtYi12YWYtbG9uZ2l0dWRpbmFsIikKcmVzdWx0c19kaXIgPC0gZmlsZS5wYXRoKGFuYWx5c2lzX2RpciwgInJlc3VsdHMiKQppbnB1dF9kaXIgPC0gZmlsZS5wYXRoKGFuYWx5c2lzX2RpciwgImlucHV0IikKZmlsZXNfZGlyIDwtIGZpbGUucGF0aChyb290X2RpciwgImFuYWx5c2VzIiwgInNhbXBsZS1kaXN0cmlidXRpb24tYW5hbHlzaXMiLCAicmVzdWx0cyIpCgojIElucHV0IGZpbGVzCnBidGFfZmlsZSA8LSBmaWxlLnBhdGgoZmlsZXNfZGlyLCAicGJ0YS50c3YiKSAjIGZpbGUgZnJvbSBhZGQtc2FtcGxlLWRpc3RyaWJ1dGlvbiBtb2R1bGUKZ2Vub21pY19wYWlyZWRfZmlsZSA8LSBmaWxlLnBhdGgoZmlsZXNfZGlyLCAiZ2Vub21pY19hc3NheXNfbWF0Y2hlZF90aW1lX3BvaW50cy50c3YiKQpwYWxldHRlX2ZpbGUgPC0gZmlsZS5wYXRoKHJvb3RfZGlyLCAiZmlndXJlcyIsICJwYWxldHRlcyIsICJvbmNvcHJpbnRfY29sb3JfcGFsZXR0ZS50c3YiKQoKIyBGaWxlIHBhdGggdG8gcGxvdCBkaXJlY3RvcnkKcGxvdHNfZGlyIDwtCiAgZmlsZS5wYXRoKGFuYWx5c2lzX2RpciwgInBsb3RzIikKaWYgKCFkaXIuZXhpc3RzKHBsb3RzX2RpcikpIHsKICBkaXIuY3JlYXRlKHBsb3RzX2RpcikKfQoKc291cmNlKHBhc3RlMChyb290X2RpciwgIi9maWd1cmVzL3NjcmlwdHMvdGhlbWUuUiIpKQpzb3VyY2UocGFzdGUwKGFuYWx5c2lzX2RpciwgIi91dGlsL2Z1bmN0aW9uLWNyZWF0ZS1iYXJwbG90LlIiKSkKYGBgCgojIFJlYWQgaW4gZGF0YSBhbmQgcHJvY2VzcwoKYGBge3IgbG9hZC1wcm9jZXNzLWlucHV0c30KcGJ0YV9kZiA8LSByZWFkcjo6cmVhZF90c3YocGJ0YV9maWxlLCBndWVzc19tYXggPSAxMDAwMDAsIHNob3dfY29sX3R5cGVzID0gRkFMU0UpICU+JSAKICBzZWxlY3QoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCwgS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRCwgY2dfbXVsdGlwbGUsIGNnX2lkLCBjZ0dGQUMsIHR1bW9yX2Rlc2NyaXB0b3IpCgpnZW5vbWljX3BhaXJlZF9kZiA8LSByZWFkcjo6cmVhZF90c3YoZ2Vub21pY19wYWlyZWRfZmlsZSwgZ3Vlc3NfbWF4ID0gMTAwMDAwLCBzaG93X2NvbF90eXBlcyA9IEZBTFNFKSAlPiUKICBsZWZ0X2pvaW4ocGJ0YV9kZiwgYnkgPSBjKCJLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lEIikpIAoKZGVzY3JpcHRvcnNfZGYgPC0gZ2Vub21pY19wYWlyZWRfZGYgJT4lCiAgZ3JvdXBfYnkoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCkgJT4lCiAgc3VtbWFyaXplKGRlc2NyaXB0b3JzID0gcGFzdGUoc29ydCh0dW1vcl9kZXNjcmlwdG9yKSwgY29sbGFwc2UgPSAiLCAiKSwpIAoKIyBWZWN0b3IgdG8gb3JkZXIgdGltZXBvaW50cwp0aW1lcG9pbnRzIDwtIGMoIkRpYWdub3NpcyIsICJQcm9ncmVzc2l2ZSIsICJSZWN1cnJlbmNlIiwgIkRlY2Vhc2VkIiwgIlNlY29uZCBNYWxpZ25hbmN5IiwgIlVuYXZhaWxhYmxlIikKCiMgQ3JlYXRlIGRmIHRvIHVzZSBmb3IgcGxvdHMKZGZfcGxvdCA8LSBnZW5vbWljX3BhaXJlZF9kZiAlPiUgCiAgbGVmdF9qb2luKGRlc2NyaXB0b3JzX2RmLCBieSA9IGMoIktpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQiLCAiZGVzY3JpcHRvcnMiKSkgJT4lIAogIG11dGF0ZSh0ZF9jZ0dGQUMgPSBjYXNlX3doZW4oZ3JlcGwoIkRlY2Vhc2VkIiwgdHVtb3JfZGVzY3JpcHRvcikgfiAieERlY2Vhc2VkIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIFRSVUUgfiB0dW1vcl9kZXNjcmlwdG9yKSwKICAgICAgICAgbG9nMTBfdG1iID0gYWJzKGxvZzEwKHRtYikpLAogICAgICAgICBjZ19pZF9raWRzID0gcGFzdGUoY2dfaWQsIEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQsIHNlcCA9ICJfIiksCiAgICAgICAgIGNnX2lkX2tpZHMgPSBzdHJfcmVwbGFjZShjZ19pZF9raWRzLCAiLyIsICJfIiksCiAgICAgICAgIGNnX2lkX2tpZHMgPSBzdHJfcmVwbGFjZShjZ19pZF9raWRzLCAiLSIsICJfIiksCiAgICAgICAgIGNnX2lkX2tpZHMgPSBzdHJfcmVwbGFjZV9hbGwoY2dfaWRfa2lkcywgIiAiLCAiXyIpLAogICAgICAgICB0dW1vcl9kZXNjcmlwdG9yID0gZmFjdG9yKHR1bW9yX2Rlc2NyaXB0b3IpLAogICAgICAgICB0dW1vcl9kZXNjcmlwdG9yID0gZmN0X3JlbGV2ZWwodHVtb3JfZGVzY3JpcHRvciwgdGltZXBvaW50cykpCgojIExldCdzIGNvdW50IG51bWJlciBvZiBzYW1wbGVzIApjb3VudF9kZiA8LSBkZl9wbG90ICU+JSAKICBncm91cF9ieSh0dW1vcl9kZXNjcmlwdG9yLCBjZ19pZCwgS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRCwgVmFyaWFudF9DbGFzc2lmaWNhdGlvbikgJT4lIAogIGRwbHlyOjpjb3VudChjZ19pZCkgCgpgYGAgCgojIERlZmluZSBwYXJhbWV0ZXJzIGZvciBwbG90cwoKYGBge3IgZGVmaW5lLXBhcmFtZXRlcnMtZm9yLXBsb3RzfQojIFJlYWQgY29sb3IgcGFsZXR0ZQpwYWxldHRlX2RmIDwtIHJlYWRyOjpyZWFkX3RzdihwYWxldHRlX2ZpbGUsIGd1ZXNzX21heCA9IDEwMDAwMCwgc2hvd19jb2xfdHlwZXMgPSBGQUxTRSkgCgojIERlZmluZSBhbmQgb3JkZXIgcGFsZXR0ZQpwYWxldHRlIDwtIHBhbGV0dGVfZGYkaGV4X2NvZGVzCm5hbWVzKHBhbGV0dGUpIDwtIHBhbGV0dGVfZGYkY29sb3JfbmFtZXMKYGBgCgojIEFsdGVyYXRpb25zIHBlciB0aW1lcG9pbnQKCmBgYHtyIHBsb3QtdGltZXBvaW50LCBmaWcud2lkdGggPSA2LCBmaWcuaGVpZ2h0ID0gNiwgZmlnLmZ1bGx3aWR0aCA9IFRSVUV9CiMgRGVmaW5lIHBhcmFtZXRlcnMgZm9yIGZ1bmN0aW9uCnhfdmFsdWUgPC0gY291bnRfZGYkdHVtb3JfZGVzY3JpcHRvcgp0aXRsZSA8LSBwYXN0ZSgiVmFyaWFudCB0eXBlcyBpbiBQQlRBIGNvaG9ydCIsIHNlcCA9ICIgIikKCiMgUnVuIGZ1bmN0aW9uCmZuYW1lIDwtIHBhc3RlMChwbG90c19kaXIsICIvIiwgIkFsdGVyYXRpb25fdHlwZV90aW1lcG9pbnRzX2JhcnBsb3RzLnBkZiIpCnByaW50KGZuYW1lKQpwIDwtIGNyZWF0ZV9zdGFja2VkX2JhcnBsb3RfdmFyaWFudChjb3VudF9kZiA9IGNvdW50X2RmLCB4ID0geF92YWx1ZSwgcGFsZXR0ZSA9IHBhbGV0dGUsIHRpdGxlID0gdGl0bGUpCnBkZihmaWxlID0gZm5hbWUsIHdpZHRoID0gNiwgaGVpZ2h0ID0gNikKcHJpbnQocCkKZGV2Lm9mZigpCmBgYAoKIyBBbHRlcmF0aW9ucyBwZXIgdGltZXBvaW50IGluIGVhY2ggY2FuY2VyIHR5cGUKCmBgYHtyIHBsb3QtY2ctaWQsIGZpZy53aWR0aCA9IDMwLCBmaWcuaGVpZ2h0ID0gMTgsIGZpZy5mdWxsd2lkdGggPSBUUlVFfQpjZ19pZF9pZCA8LSBhcy5jaGFyYWN0ZXIodW5pcXVlKGNvdW50X2RmJGNnX2lkKSkgCmNnX2lkX2lkIDwtIHNvcnQoY2dfaWRfaWQsIGRlY3JlYXNpbmcgPSBGQUxTRSkKY2dfaWRfaWQKCiMgRGVmaW5lIHBhcmFtZXRlcnMgZm9yIGZ1bmN0aW9uCnhfdmFsdWUgPC0gY291bnRfZGYkS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRAp0aXRsZSA8LSBwYXN0ZSgiVmFyaWFudCB0eXBlcyBpbiBQQlRBIGNvaG9ydCBhY3Jvc3MgY2FuY2VyIGdyb3VwcyIsIHNlcCA9ICIgIikKCiMgUnVuIGZ1bmN0aW9uCmZuYW1lIDwtIHBhc3RlMChwbG90c19kaXIsICIvIiwgIkFsdGVyYXRpb25fdHlwZV90aW1lcG9pbnRzX2NnX2lkX2JhcnBsb3RzLnBkZiIpCnByaW50KGZuYW1lKQpwIDwtIGNyZWF0ZV9zdGFja2VkX2JhcnBsb3RfdmFyaWFudF9jZ19pZChjb3VudF9kZiA9IGNvdW50X2RmLCB4ID0geF92YWx1ZSwgcGFsZXR0ZSA9IHBhbGV0dGUsIHRpdGxlID0gdGl0bGUpCnBkZihmaWxlID0gZm5hbWUsIHdpZHRoID0gMzAsIGhlaWdodCA9IDE4KQpwcmludChwKQpkZXYub2ZmKCkKCmBgYAoKIyBBbHRlcmF0aW9ucyBwZXIgdGltZXBvaW50IGluIGVhY2ggY2FuY2VyIHR5cGUgZGVmaW5lZCBieSBjZ0dGQUMKCmBgYHtyIHBsb3QtY2dHRkFDLW4taW5kaXZpZHVhbC1wbG90cywgZmlnLndpZHRoID0gNSwgZmlnLmhlaWdodCA9IDUsIGZpZy5mdWxsd2lkdGggPSBUUlVFfQpkZl9wbG90X2NnR0ZBQyA8LSBkZl9wbG90ICU+JQogIGZpbHRlcighaXMubmEodGltZXBvaW50c19tb2RlbHMpKSAlPiUgCiAgYXJyYW5nZSh0aW1lcG9pbnRfY2dHRkFDX24pICU+JSAKICBncm91cF9ieSh0dW1vcl9kZXNjcmlwdG9yLCBjZ0dGQUMsIHRpbWVwb2ludF9jZ0dGQUNfbiwgS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRCwgVmFyaWFudF9DbGFzc2lmaWNhdGlvbikgJT4lIAogIGRwbHlyOjpjb3VudCh0aW1lcG9pbnRfY2dHRkFDX24pCgpjZ0dGQUNfaWQgPC0gYXMuY2hhcmFjdGVyKHVuaXF1ZShkZl9wbG90X2NnR0ZBQyRjZ0dGQUMpKSAKY2dHRkFDX2lkIDwtIHNvcnQoY2dHRkFDX2lkLCBkZWNyZWFzaW5nID0gRkFMU0UpCmNnR0ZBQ19pZAoKIyBEZWZpbmUgcGFyYW1ldGVycyBmb3IgZnVuY3Rpb24KeF92YWx1ZSA8LSBkZl9wbG90X2NnR0ZBQyRLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lECnRpdGxlIDwtIHBhc3RlKCJWYXJpYW50IHR5cGVzIGluIFBCVEEgY29ob3J0IGFjcm9zcyBjZ0dGQUMiLCBzZXAgPSAiICIpCgojIFJ1biBmdW5jdGlvbgpwIDwtIGNyZWF0ZV9zdGFja2VkX2JhcnBsb3RfdmFyaWFudF9jZ0dGQUMoY291bnRfZGYgPSBkZl9wbG90X2NnR0ZBQywgeCA9IHhfdmFsdWUsIHBhbGV0dGUgPSBwYWxldHRlLCB0aXRsZSA9IHRpdGxlKQoKYGBgCgojIEFsdGVyYXRpb25zIHBlciB0aW1lcG9pbnQgaW4gZWFjaCBjYW5jZXIgdHlwZSBhbmQgdGltZXBvaW50IG1vZGVsCgpgYGB7ciBwbG90LXRpbWVwb2ludC1tb2RlbCwgZmlnLndpZHRoID0gMTUsIGZpZy5oZWlnaHQgPSAxMiwgZmlnLmZ1bGx3aWR0aCA9IFRSVUV9CnRtX2RmX3Bsb3QgPC0gZGZfcGxvdCAlPiUKICBmaWx0ZXIoIWlzLm5hKHRpbWVwb2ludHNfbW9kZWxzKSkgJT4lIAogIGdyb3VwX2J5KHR1bW9yX2Rlc2NyaXB0b3IsIGNnX2lkLCB0aW1lcG9pbnRzX21vZGVscywgS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRCwgVmFyaWFudF9DbGFzc2lmaWNhdGlvbikgJT4lIAogIGRwbHlyOjpjb3VudCh0aW1lcG9pbnRfY2dHRkFDX24pCgp0bSA8LSBhcy5jaGFyYWN0ZXIodW5pcXVlKHRtX2RmX3Bsb3QkdGltZXBvaW50c19tb2RlbHMpKQp0bSA8LSBzb3J0KHRtLCBkZWNyZWFzaW5nID0gRkFMU0UpCnRtCgojIExvb3AgdGhyb3VnaCB2YXJpYWJsZQpmb3IgKGkgaW4gc2VxX2Fsb25nKHRtKSl7CiAgcHJpbnQoaSkKICBkZl9zdWIgPC0gdG1fZGZfcGxvdCAlPiUKICAgICAgZmlsdGVyKHRpbWVwb2ludHNfbW9kZWxzID09IHRtW2ldKQogIAogICAjIERlZmluZSBwYXJhbWV0ZXJzIGZvciBmdW5jdGlvbgogIHhfdmFsdWUgPC0gZGZfc3ViJEtpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSUQKICB0aXRsZSA8LSBwYXN0ZSh0bVtpXSkKICAKICAjIFJ1biBmdW5jdGlvbgogIHAgPC0gY3JlYXRlX3N0YWNrZWRfYmFycGxvdF92YXJpYW50X2NnX2lkKGNvdW50X2RmID0gZGZfc3ViLCB4ID0geF92YWx1ZSwgcGFsZXR0ZSA9IHBhbGV0dGUsIHRpdGxlID0gdGl0bGUpCiAgCn0KYGBgCgoKYGBge3IgZWNobz1UUlVFfQpzZXNzaW9uSW5mbygpCmBgYAo=